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PsJ 1. Introduction 

CN Sub Shot Noise (SSN) imaging of weak object by exploiting Parametric Down Con- 

verted light, i.e. to obtain the image of a weak absorbing object with a level of noise 
below the minimum threshold that is unavoidable in the classical framework of light 
detection, represents a very interesting technological development |i;H'"92|dit|and[ ^ 
precise characterization of PDC speckle structure in dependence of pump beam 
parameters is a fundamental tool for this application. Indeed, it is fundamental to 
5^ set the dimension of the modes coherence areas with respect to the dimension of 

the pixels of the ccd camera used for acquiring the images. 

Albeit a theory of speckles structure in PDC l '"9 | ™°* l has been developed and 
some results in a specific regime have been collected a systematic study of this 
structure is still missing. In this paper we present a first set of data addressed to 
this purpose. 
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2. Theory 

The process of SPDC presents a large bandwidth in the spatial frequency domain, 
that is particularly useful when studying spatial quantum correlations 1^. Any pair 
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of transverse modes of the radiation, characterized by two opposite transverse mo- 
menta q and — q, are correlated in the photon number, i.e. they contain, in an ideal 
situation, the same number of photons. In the far field zone, the single transverse 
mode is characterized by a coherence area, namely the uncertainty on the emission 
angle (tan?? = Ag/27r, A being the wavelength) of the twin photons. It derives 
from two effects that participate in the relaxation of the phase matching condition. 
On the one side the finite transverse dimension of the gain area, coinciding with 
the pump radius Wp at low parametric gain. On the other side the finite longitu- 
dinal dimension of the system, i.e. along the pump propagation direction, that is 
generally given by the crystal length I. If the first dominates, the coherence area is 
related to the Fourier transform of the pump transverse profile, i.e. Xcoh oc 1/wp. 
If the second dominates, the coherence area is of the order of {I ■ ?9)~^ for small 
emission angle. The appearance of the emission is a speckled structure in which the 
speckles have, roughly, the dimension of the coherent area and for any speckle at 
position q there exists a symmetrical one in — q with equal intensity (see fig. |2|. 

Omitting some unessential constants, the Hamiltonian describing the three fields 
parametric interaction is 

Hj{t) cc f x'^' E(+)(r,i) Ei+\r,t) E(-)(r,i) d\ + h.c (1) 
Jv 

The pump depletion due to the down-conversion and the absorbtion is indeed of 
small entity, unless extremely high intensity laser sources are used. We shall there- 
fore work in the parametric approximation, that treats the pump as a classical 
monochromatic field propagating linearly along a certain z direction inside the 
crystal and having an amplitude transverse profile Ap(p), i.e. 

Ep{v,t)ocApip)e-'^''-'-^-*\ (2) 

where p is the coordinate vector in the transverse plane to the propagation direction 
z. 

The down-converted fields 1 and 2 are quantized. Their positive- and negative- 
frequency part E'^+^(r,t) and E^^^(r,t) is given as an expansion in plane-wave 
modes and we find convenient to express them separating the sum over the wave- 
vector into the sum over its transverse component q and the frequency uj. Thus, we 
have 

e(+) (r, i) oc ^ e^(k..^-^a) g^q.•P - (q^^ (3) 

where i = 1,2. The third component kiz of the i-th field wave vector is expressed 
in terms of the q^ and oji because of the relations 

hz = \J kf - g2 and h = (4) 

Introducing Eq.s Q and ^ in the Hamiltonian (jlJ we have 
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e*^"''-"^-"^)* ai(qi,c^i)a2(q2,c^2) + h.c. (5) 



qi,i^i q2,i^2 



Here I is the length of the crystal, while x ly is the area of its transverse surface. 

The integral in dz gives a contribution proportional to I ■ sine [(Afc /) /2)] where 
Afc = fciz + k2z — kp. The double integral on the transverse surface of the crystal 
gives the Fourier transform of the pump transverse profile if the crystal is large 
compared to it. Supposing a gaussian pump Ap(p) = ApC^'' the Hamiltonian 
([5]) becomes 

"Afc(qi,q2, wi,u;2) ■ 



Hi oc 



E 



qi,i»Ji q2,i^2 



g ■ smc 



-(qi+q2)-' 



g^K^^i-^2)t2^(q^^^) S2(-q,t^2) + h.c. 



(6) 



where we have introduced the dimensionless factor g (x x*"^"* ■ I ■ Ap, usually 
referred to as parametric gain. Its value determines the number of photons that 
are generated in the down conversion process in mode pairs that are well-phase 
matched. The evolution of the quantum system guided by Hamiltonian ([6|, in the 
case of relatively high gain regime, requires a numerical solution and it is discussed 
in detail inE^ Anyway, in the first order of the perturbation theory {g <^ 1), the 
quantum state of the scattered light has the entangled form 



IV") — |vac) + exp 



Hidt 



|0) 



|vac)+ ^ ^F(qi,q2,f^)|lq„o)|lq„- 



(7) 



^(qi,q2,f^) = 9 ■ sine 



Ak{qi,q2,n) -l 



UJi — LUp/2 + ri, U>2 — Wp/2 — fl. 



(8) 



The coherence area, in the limit of low parametric gain g, can be estimated 
by the angular structure of the coincidence probability |-F(qi,q2)| at some fixed 
frequency f2. As mentioned before, now is clear that we deal with two functions 
that enter in the shaping of the coherence area: the sine function and the Fourier 
transformed of the gaussian pump profile. Since they are multiplied, the narrower 
determines the dimension of the area. By expanding linearly the longitudinal wave 
detuning around the exact matching point Afc(qo,qo,fi) according to relations 
Q, the sine function turns out to have a Half Width Half Maximum of Aq = 
2, 78 /{I tan 1}). The HWHM of the gaussian function, appearing in ([s]), is 5g = 
y/2 ln(2) /wp. Concerning our experiment, we consider so small emission angles i? 
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and large enough pump radius Wp, that we always work in the region Sq/Aq < 1. 
Therefore, in principle, the dimension of the coherence area is only determined by 
the pump waist. 

When moving to higher gain regime, the number of photon pairs generated in 
the single mode increases exponentially as cx sinh^( g) i.e. a large number of photons 
is emitted in the coherence time along the direction i?. In this case, also the pump 
amplitude becomes important in the determination of the speckles dimension. As 
described in this can be explained by a qualitative argumentation: inside the 
crystal, the cascading effect that causes the exponential growth of the number of 
generated photons is enhanced in the region where the pump field takes its highest 
value, i.e. close to the center of the beam. Thus, in high gain regime, most of the 
photon pairs are produced where the pump field is closed to its peak value. As 
a result the effective region of amplification inside the crystal becomes narrower 
than the beam profile. Thus, in the far field one should consider the speckles as the 
Fourier transform of the effective gain profile, that being narrower, produces larger 
speckles. 

A further fundamental consideration for the practical implementation is that in 
high gain regime, instead of measuring the coincidences between two photons by 
means of two single photo-detectors, one collects a large portion of the emission 
by using for instance a CCD array with a certain fixed exposure time. Within 
this time several photons are collected by the single pixels and the result is an 
intensity pattern, having the spatial resolution of the pixel. The coherence area can 
be evaluated by the cross-correlation between the signal's and the idler's intensity 
patterns. We can define also the auto-correlation function of the signal intensity 
pattern itself, since in the single transverse mode of the signal arm there are many 
photons. To be precise the speckle's dimension is better related to the spread of 
this function, although the two functions, the cross- and the auto- correlation, 
present the same behaviour with respect to the pump parameters, see liMi From 
the experimental view-point it is convenient to study the auto-correlation because 
of the higher visibility that allows a more accurate estimation of its size. 

3. Experiment 

In our setup, Figjl] a type II BBO non-linear crystal = 1 cm) is pumped by the 
third harmonic (wavelength of 355 nm) of a Q-switched Nd:Yag laser. The pulses 
have a duration of 5ns with a repetition rate of 10 Hz and a maximum energy, at 
the selected wavelength, of about 200 mJ. The pump beam crosses a spatial filter 
(a lens with /=50 cm and an iris of 250 /im of diameter), in order to eliminate the 
non-gaussian components and to coUimate it before the crystal. The diameter of the 
pump beam entering the crystal is varied, when necessary, by changing the distance 
between two lenses (a biconvex and a biconcave) placed after the spatial filter. After 
the crystal, the pump is stopped by a UV mirror, transparent to the visible, and 
by a low frequency-pass filter. The down converted photons (signal and idler) pass 
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Fig. 1. Experimental setup. A triplicated Nd-Yag laser beam, after spatial filtering, produces 
type II PDC in a BBO crystal, which is then measured, after an interference filter and pump 
elimination, by a CCD camera. 



through a lens of 5 cm of diameter (/ = 10 cm) and an interference fiher centered 
at the degeneracy A=710 nm (lOnm bandwidth) and finally measured by a CCD 
camera. We used a 1340X400 CCD array, Princeton Pixis:400BR (pixel size of 20 
/Ltm), with high quantum efficiency (80%) and low noise (5 electrons/pixel). The far 
field is observed at the focal plane of the lens in a / — / optical configuration, that 
ensures that we image the Fourier transform of the crystal exit surface. Therefore 
a single transverse wavevector q is associated to a single point x = (A//27r)q in 
the detection plane. The CCD acquisition time is set to 90 ms, so that each frame 
corresponds to the PDC generated by a single shot of the laser. 

Looking at the images, we can appreciate the speckled structure and a certain 
level of correlation of the speckles intensity between the signal and idler arms (Fig. 
[2]). Let us define iV/j(x) the intensity level, proportional to the number of photons, 
registered by the pixel in the position x of the region R. SNj-i{x.) = Nfj{x.) — 
{Nii{x)) is the fluctuation around the mean value that is estimated as (Affl(x)) = 
(1/n) A^fl(x), with n the number of pixels. We evaluate the normalized spatial 
cross-correlations of the intensity fluctuations in an arbitrary region Ri , belonging 
to the signal portion of the image, and in the symmetric region i?2 (see Fig. [2| 
belonging to the idler portion: 
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where ^ is the displacement vector that assumes discrete values. C12 reaches 
a peak of about 0,9 in Fig. [2]-a, that indicates a good level of spatial correlation 
between signal and idler. In figure [2]-b the value is around 0,6 because here we did 
not put the interference filter in front of the CCD camera, allowing more background 
light to enter, and the signal and idler component are not separated. It is worth to 
emphasize that Ci2(0 "^o^ ^ index of the correlation at the quantum level and 
it can not be used to discriminate the SSN condition. We are interested mainly 
in the width of the peak, that indicate the coherence area. On the other hand, 
the SSN condition is checked by a measurement of ct^ =< {Nn-^ — Nj^^)'^ > — < 
— Nfi^ >^< 1: our preliminary data show that this result is reached for some 
images without background subtraction; a detailed discussion of this investigation 
is postponed to a forthcoming paper. 

Since the experimental cross-correlation is characterized by a non-optimal vis- 
ibility, due to the losses in optical paths and noise, we preferred to estimate the 
speckle's size, by evaluating the auto-correlation of a single region R of the signal 
intensity pattern. 

First off all we investigate the gain region in which we are working, in order 
to ensure the possibility to reach a sufficient non-linear cascading effect in the 
photon pairs production inside the crystal. Fig. [3] shows the mean photon number 
{Nji(x)) as function of the pump power PW. Any point is averaged over several 
tenths of frames in order to reduce the uncertainty. Our laser presents in fact 20% 
fluctuations of the power from pulse-to-pulse. Since the mean number of photons is 
proportional to smh^^g), the fluctuations of the pump generate large fluctuations 
in the photons number from framc-to-framc. The diameter of the pump is fixed 
and the power is varied by the delay between the Q-switch turn-on and the lamp 
flash. The calibration curve delay-power has been measured by a power meter and 
we observed a reproducibility with uncertainty around 10%. 

In Fig. |3]-a the pump diameter is around 1,3mm and in Fig. [3]-b is 0,95mm. 
Although the pump power used for the case (b) is smaller than case (a), the intensity 
results higher in the case (b) because of the reduced radial dimension. It must be 
noticed that we are constrained in the range of intensities of the pump. For high 
intensities we are limited by the damage threshold of optical components, for low 
intensities by the visibility of the speckle structure because we collect a lot of 
temporal modes in the same frames. 

The data are fitted by the equation {Nfj) ^ k ■ smh'^{ay/p). k depends from 
the number of modes while a\/P = g. The experimental values, mediated on three 
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Fig. 2. Spatial correlation between signal and idler down converted light with interference filter 
around degeneracy{a) and without a narrow filter (b). In each raw the first inset is a ccd image, 
the second a plot of the cross correlation function. The axis report the pixels position in the CCD 
array. Ci2(^). 



different acquisitions, are k — 31.48 and a = 1.91 for the case (a). In the case (b) 
we obtain k = 1, 10 and a — 4,87; the higher value of cr is due to the increased 
intensity. Therefore, in the case (a) we have a gain g from 1,5 to 3,5 and in the case 
(b) from 1,9 to 5. Thus, we are in a non- linear regime where we should expect a 
dependence of the speckles size from the amplitude of the pump. 

Fig. |3]-c and|3]-d show the trend of the radius of the coherence area with the pump 
power, again at the two different diameters of the pump. Actually we observe an 
increasing of the radius, that is predicted by the considerations exposed in section 
[2] We consider, just as indicative, a linear fit y = a{x — b), obtaining a = 1, 25 and 
b — —0,63 for (c) and a — 3.7 and 6 = —0,27 for (d). The higher slope of (c) is 
qualitatively explained by the higher gain. 

Finally, we investigate the dependence of the radius of the speckles from the 
pump diameter, as shown in figure |4] We fix the power of the laser to 0,78 MW and 
we change the diameter varying the distance between the two coUimating lenses 
(see figure PI . The two curves differ just in the estimation of the pump diameter: 



September 30, 2008 9:6 WSPC/INSTRUCTION FILE specklel 



{t>> 



Fig. 3. The top of the figure represents the mean number of photon/pixel collected by the CCD 
in one shot of the laser (5ns), function of the pump pulse power in MW. The bottom represents 
the dependence of the speckle radius (in pixel) from the pump power. In (a) and (c) the pump 
diameter is 1,3mm, while in (b) and (d) is 0.95mm. 



in one case we measured it with an impact paper (IP) and in the other case with 
a CCD. We perform a fit in the form y ^ a ■ obtaining a = (8, 1 ± 0, 1) and 
b = (-3, 61 ± 0, 09) for the IP curve and a = (3, 22 ± 0, 07) and b = (-3, 73 ± 0, 09) 
for the CCD curve. Despite the different way of pump size estimation, reflected in 
the value of the parameter a, the coefficients b are compatible. The theory, in low 
gain regime, provides that the radius of the speckles is proportional to the inverse of 
the pump size (wp), i.e. b — —1. The estimated value of b, in our case, confirms the 
role of the high gain regime in the speckle size. In fact, together with the reduction 
of the pump diameter, the gain increases, and thus the effective gain area is more 
reduced again. This effect impresses upon the speckles size a stronger dependence 
with respect to the pump size. 
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Fig. 4. Observed dependence of the radius (in pixels) of the speckles from the pump diameter. 
The two curves differs in the estimation of the pump diameter: in one case it was measured with 
an impact paper (IP) and in the other case with a CCD. 



4. Conclusion 

In this paper we have presented a first experimental study of the size of the co- 
herence area in PDC in the high gain regime. We have shown that the speckles 
present, not only a dependence on the pump diameter, as in the usual low gain 
regime, but also a strong dependence from the pump intensity. This result, i.e. the 
understanding of the behaviour of the coherence area in the high gain regime, is 
fundamental for the innovative application in the field of quantum imaging but 
also in general for quantum metrology, quantum information, etc. A comparison 
with these data will represent a challenge for the theoretical models describing this 
phenomenon. 
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